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Abstract.  Infrared  radiative  excitation  in  non-local  thermodynamic  equilibrium  (non-LTE) 
regions  of  the  Earth's  atmosphere  for  the  V3-mode  vibrationally  excited  states  of  CO2  under 
sunlit  conditions  and  the  resulting  4.3-pm  limb  radiance  are  calculated  using  a  line-by-line  (LBL) 
radiative  transfer  model.  Excited-state  population  densities  and  the  corresponding  vibrational 
temperature  profiles  are  calculated  for  the  important  emitting  states  using  a  model  which  includes 
radiative  absorption  and  emission  as  well  as  various  collisional  processes.  The  quenching  of 
0 {lD)  by  N2  has  a  greater  impact  on  these  population  densities  than  has  been  previously 
reported  in  the  literature.  Integrated  radiance  in  a  limb  view  for  the  4.3-pm  bands  is  calculated 
trom  the  model  and  compared  with  sunlit  earthlimb  measurements  obtained  by  the  Spectral 
Infrared  Rocket  Experiment  (SPIRE).  Solar  pumping  is  the  dominant  excitation  process  for  the 
4.3-pm  emitting  states  in  the  daytime.  The  major  contribution  to  the  total  limb  radiance  for 
tangent  heights  ot  55-95  km  is  made  by  the  fluorescent  states  at  approximately  3600  cm** 
which  absorb  sunlight  at  2.7  pm  and  then  emit  preferentially  at  4.3  pm .  The  predicted  radiance 
^isrin  good  agreement  with  the  SPIRE  measurements  for  all  tangent  heights  in  the  50-  to  130-km 
>ange.  This  is  the  first  detailed  comparison  of  results  of  a  full  line-by-line  non-LTE  radiative 
transfer  calculation  with  4.3-pm  earthlimb  radiance  data. 

1.  Introduction 

In  recent  years,  there  has  been  increasing  interest  in 
radiative  transfer  in  the  infrared  spectral  region  under  non¬ 
local  thermodynamic  equilibrium  (non-LTE)  conditions  in  the 
Earth's  atmosphere  [e.g..  Kumer .  1977;  Shved  et  al,  1978; 

Kumer  and  James,  1982;  Gordiets  et  al. ,  1982;  Dickinson , 

1984;  Lopez-Puertas  et  al,,  1986a. b;  Lopez-Puertas  and 
Taylor,  1989;  Kittepov  et  al.,  1991;  Wintersteiner  et  al., 

1992].  These  non-equilibrium  conditions  generally  exist  in 
the  Earth's  mesosphere  and  thermosphere,  where  the  time 
between  collisions  is  not  much  less  than  the  radiative 
lifetime.  Non-LTE  infrared  emissions  play  an  important  role 
in  establishing  the  global  heat  balance,  structure,  and 
dynamical  properties  of  the  middle  atmosphere. 

In  the  present  paper.  4.3  pm  infrared  radiation  associated 
with  the  v3  (asymmetric  stretch)  vibrational  mode  of  C02  is 
considered  in  the  non-LTE  region  of  the  quiescent  atmosphere 
under  day  lit  conditions.  A  future  paper  will  deal  with  the 
nighttime  emission.  A  set  of  codes  known  collectively  as 
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atmospheric  radiance  code  (ARC)  is  used  to  calculate  (1)  the 
vibrational  populations  of  states  that  emit  at  4.3  pm  and  (2) 
the  resulting  earthlimb  radiance.  The  model  calculations  are 
compared  to  the  daytime  earthlimb  radiance  data  from  the 
Spectral  Infrared  Rocket  Experiment  (SPIRE)  [Stair  et  aL, 
1985]  and  it  is  found  that  the  model  results  are  in  excellent 
agreement  with  the  measurements. 

One  noteworthy  aspect  of  the  ARC  codes  is  that  all  infrared 
radiative  processes  are  treated  with  line-by-line  (LBL) 
algorithms.  In  particular,  the  vibrational  populations,  largely 
determined  by  the  transfer  of  radiation  within  the  atmosphere 
and  by  the  rate  of  absorption  of  solar  flux,  are  calculated  with 
the  RAD  code  [Wintersteiner  et  al.,  1992].  The  solar  flux 
absorption  rates  come  from  a  separate  LBL  code,  SABS 
[Wintersteiner  and  Joseph .  1986],  and  the  earthlimb  radiance 
is  subsequently  calculated  from  the  level  populations  with  a 
third  code,  NLTE  [Wintersteiner  and  Sharma.  1985].  This 
approach  has  previously  been  used  to  model  the  CCb  15-pm 
emission  seen  by  the  SPIRE  experiment  [ Wintersteiner  et  al., 
1992]. 

Comparisons  between  model  calculations  and  4.3  pm 
earthlimb  radiance  measurements  have  been  carried  out 
previously  by  Kumer  and  James  [1982]  using  the  SPIRE  data 
set  and  by  Lopez-Puertas  and  Taylor  [1989]  using 
Stratospheric  and  Mesospheric  Sounder  (SAMS)  data,  but 
neither  of  these  studies  calculated  the  radiative  transfer  using  a 
full  line-by-line  treatment.  Kumer  and  James  [1982]  used  a 
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band  model  to  calculate  the  radiative  transfer.  Lopei-Puertas 
and  Taylor  [1989]  used  line-by-line  frequency  integration  to 
calculate  the  transmittances  of  the  pressure-modulated  cell  used 
in  the  SAMS  experiment,  but  they  used  a  narrowband  model  to 
calculate  the  radiative  transfer.  In  the  present  calculation  the 
radiative  transfer  is  calculated  by  integrating  over  frequency 
for  each  spectral  line  of  a  given  rovibrational  band,  using 
line-specific  parameters  taken  from  the  HITRAN  molecular 
absorption  parameters  database  [Rothman  et  ai ,  1987]  for 
each  line.  The  altitude  dependence  of  the  Voigt  emission  and 
absorption  line  shapes  as  well  as  the  altitude  dependence  of 
the  line  strengths  are  taken  into  account  for  each  line  in  each 
atmospheric  layer.  Thus  the  exchange  of  photons  among 
layers  is  calculated  very  accurately. 

Section  2  contains  an  outline  of  the  model,  including  the 
coupling  scheme  used  for  the  relevant  vibrational  states,  a 
discussion  of  the  model  atmosphere  used  for  the  SPIRE 
simulation,  and  a  description  of  the  SPIRE  data  set.  Section 
3.1  contains  the  results  of  the  calculations  for  the  vibrational 
temperatures  of  all  the  states  considered.  In  section  3.2  the 
results  of  the  limb  radiance  calculations  are  presented  and 
compared  with  the  sunlit  portion  of  the  SPIRE  data  set. 
Section  4  contains  a  comparison  of  the  results  with  those  of 
earlier  studies,  as  well  as  a  discussion  of  approximations  and 
uncertainties.  In  section  5  the  main  points  are  summarized. 

2.  Outline  of  the  Model 

2.1.  Coupling  Scheme 

The  vibrational  energy  levels  of  C02  used  in  the  present  set 
of  calculations  are  shown  in  Figure  1.  The  notation  for  the 
levels  is  that  of  the  HITRAN  molecular  absorption  parameters 
database  [Rothman  et  a!.,  1987,  1992],  as  explained  by 
McClatchey  et  al.  [1973].  The  V3  transitions  near  4.3  pm 


involving  emission  of  one  v3  quantum  are  very  strong,  with 
Einstein  A  coefficients  -400  s'1,  rather  than  the  -10  s'1 
values  typical  for  allowed  infrared  rovibrational  transitions. 
The  lowest-lying  4.3-pm  transition  is  the  v3  fundamental, 
00011-00001.  The  lowest-lying  hot  band,  01111-01101,  is 
also  important  because  of  the  significant  probability  that  the 
lower  state  is  occupied.  Other  important  4.3-pm  bands 
originate  in  higher-lying  states  having  Vj  and/or  v2 
excitation  as  well  as  v3  excitation.  The  states  in  groups  1  and 
2  (see  Figure  1)  are  populated  by  transitions  at  4.3  pm  and  2.7 
pm,  while  those  in  group  3  are  populated  by  2.0-pm 
transitions  from  the  ground  state.  In  other  words,  absorption 
of  solar  flux  at  2.7  and  2.0  pm  is  followed,  with  high 
probability,  by  emission  at  4.3  pm,  and  any  simulation  of 
daytime  conditions  must  properly  account  for  these  solar 
fluorescent  processes  [James  and  Kumer,  1973;  Sharma  and 
Wintersteiner,  1985]. 

It  is  important  to  consider  the  contributions  of  three  of  the 
minor  isotope  bands  in  addition  to  the  bands  of  the  principal 
l2C,602  isotope.  The  convention  used  here,  following 
McClatchey  et  ai  [1973],  is  to  designate  the  major  isotope  by 
the  digits  626  and  the  important  minor  isotopes  by  636,  628, 
and  627  to  represent  substitutions  of  13C,  180.  and  170, 
respectively.  Table  1  lists  ail  of  the  bands  considered  for  this 
study  along  with  some  of  their  important  properties.  Line 
parameters  for  each  of  the  rovibrational  transitions  have  been 
taken  from  the  1986  version  of  the  HITRAN  database 
[Rothman  et  ai ,  1987], 

A  coupling  mechanism  that  connects  all  C02  states  having 
v3  excitations  is  the  near-resonant  V-V  process  whereby  N2 
vibration  is  exchanged  with  quanta  of  asymmetric  stretch 
according  to 

C02  (v  1  v2v3)  +  N2  (v  =  0)  <=> 

C02(V|V2v3-l)  +  N2(v=  1),  (1) 


Energy 
(cm' l) 


Vibrational  Modes 


Figure  1.  C02  vibrational  levels  and  4.3-pm  transitions  considered  in  this  study.  Absorption  at  2.7  and 

2.0  pm  is  not  shown. 
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Table  1.  CO2  Bands  Included  in  the  Model  for  Emission  at  4.3  \im 


Isotopea 

Upper 

Stateb 

Lower 

State*3 

Band 

Origin 

cm"* 

Transition 

Wavelength 

fim 

Einstein  A 
Coefficient 
s*  1 

626 

00011 

00001 

2349.14 

4.26 

436.0 

01111 

01101 

2336.63 

Group  1 

4.28 

426.3 

10012 

10002 

2327.43 

4.30 

415.4 

02211 

02201 

2324.14 

4.30 

418.5 

10011 

10001 

2326.60 

Group  2 

4.30 

412.7 

11112 

11102 

2315.24 

4.32 

408.1 

03311 

03301 

2311.67 

4.33 

412.5 

11111 

11101 

2313.77 

Group  3 

4.32 

404.1 

20013 

20003 

2305.26 

4.34 

398.8 

20012 

20002 

2306.69 

4.34 

389.6 

20011 

20001 

2302.52 

4.34 

391.6 

04411 

04401 

2299.21 

4.35 

402.5 

12212 

12202 

2302.96 

4.34 

400.9 

12211 

12201 

2301.05 

4.35 

395.5 

636 

00011 

00001 

2283.49 

4.38 

371.5 

01111 

01101 

2271.76 

Group  I 

4.40 

368.1 

10012 

10002 

2261.91 

4.42 

364.8 

02211 

02201 

2260.05 

4.42 

361.7 

10011 

10001 

2262.85 

4.42 

367.2 

628 

00011 

00001 

2332.11 

4.29 

388.2 

627 

00011 

00001 

2340.01 

4.27 

395.7 

a  626  means  l60,2C!6C  and.  similarly,  for  636,  628,  and  627. 
k  Listed  in  Benedict  (HITRAN)  notation  [McClatchey  et  al. ,  1973]. 


where  any  isotope  of  C02  may  be  involved.  The  rate  constant 
for  this  process  was  originally  calculated  by  Sharma  and  Brau 
[1969].  It  has  been  measured  for  the  626  isotope  [Rosser  et 
al 1969;  Gueguen  et  al..  1975;  Inoue  and  Tsuchiva .  1975], 
and  the  same  value  is  assumed  for  the  other  isotopes.  Its  value 
of  5-6  x  10” 13  cm3/mol  sec  at  mesospheric  and  lower 
thermospheric  temperatures  is  great  enough  that  the  process  is 
a  determining  factor  for  the  vibrational  populations  of  both 
COi  and  N2,  which  therefore  cannot  be  calculated 
independently  [Kumer  and  James .  1974].  Thus  as  the  first  step 
in  these  calculations,  the  00011  states  for  the  four  isotopes 
626,  636,  628,  and  627  are  coupled  together  along  with  the 
N2(v=1)  state.  Population  densities  and  vibrational 
temperatures  for  these  five  states  are  calculated  simultaneously 
using  the  line-by-line  technique  of  the  RAD  code,  considering 
as  production  and  loss  processes  absorption  of  sunlight  and  of 
earthshine  (that  is.  photons  originating  elsewhere  in  the 
atmosphere),  spontaneous  emission,  and  the  coliisional 
processes  l- 10  indicated  in  Table  2.  All  absorbed  and  emitted 
photons  are  in  the  4.3-pm  band.  Populations  of  the  group  1 
states  are  taken  as  known  quantities  for  this  calculation:  they 
are  used  for  the  determination  of  the  N2(v=l)  vibrational 
temperature  because  this  state  is  coupled  to  the  group  1  states 
as  well  as.  to  the  000 IT  states;  Coefficients  for  earthshine 
absorption  are  calculated  by  the  RAD  algorithm.  For  the 


daytime  conditions  considered  here,  direct  solar  pumping  at 
4.3  Jim  is  accounted  for  by  use  of  solar  flux  absorption 
coefficients  from  the  line-by-line  SABS  code  [Wintersteiner 
and  Joseph ,  1986].  Rate  constants  for  the  coliisional 
interactions  are  the  same  as  those  used  by  Lopez-Puertas  et  al. 
[1986b].  The  quenching  of  0([D)  by  N2  (reaction  (7)  of  Table 
2)  is  assumed  to  take  place  with  25%  efficiency  [Harris  and 
Adams ,  1983]  which  means  that  on  average  25%  of  the  energy 
available  from  the  reaction  appears  as  N2  vibrational  energy. 
A  single  collision  can  produce  N2  with  as  many  as  seven 
vibrational  quanta:  we  assume  that  multiquantum  states  quickly 
relax  to  an  equivalent  number  of  singly  excited  N2  molecules. 

The  non-LTE  populations  of  states  of  the  C02(V|V2)  (bend- 
stretch)  manifold  which  are  used  in  these  calculations  are  also 
regarded  as  known  quantities  and  are  assumed  to  be 
independent  of  the  state  of  excitation  of  states  with  V3  quanta 
which  radiate  at  4.3  (im.  They  are  calculated  by  procedures 
described  by  Wintersteiner  et  al.  [1992],  and  this  assumption 
is  discussed  therein.  The  states  of  the  VjV2  manifold  include 
(see  Figure  1)  the  01101  state  having  one  v2  quantum,  2v2 
states  with  energy  near  1335  cm'1  (02201,  10001,  and 
10002),  and  3v2  states  which  lie  near  2000  cm'1  (03301, 
11101,  and  11102).  It  is  assumed  that  the  2v2  states  are  in 
mutual  equilibrium,  with  the' same  assumption  for  the  3v2 
states.  The  vibrational  temperature  profiles  used  here  for  the 
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Table  2.  Reactions  Included  in  the  Model  for  the  'C02(0001  1),  N2(l),  and  'CO2(01 1 1 1)  Populations 


Reaction 

Number 

Reaction3 

Rate  Constant13 
cm3/mol  s 

1 

'C02(0001 1)  +  N2(0)  »  'C02(00001)  +  N2(l)  +  19.2  cm*1 

5.0  x  10-13(300/T),/2 

2 

'C02(00011)  +  M  <=>  'C02(03301)  +  M +  345.9  cm'1 

2.22  x  10‘15 

+  1.22  x  10"10  exp(-76.75/T°-33) 

3 

'C02(00011)  +  O2(0)  'CO2(01 101)  +  02(1)  +  125.4  cm"1 

3  x  lO'15  (1  +  0.02(T  -  210)) 

4 

'C02(00011)  +  0  »  'C02(03301)  +  O  +  345.9  cm’1 

2  x  10‘13  (T/300) !/2 

5 

N2(l)  +  O2(0)  <=>  N2(0)  +  O2(l)  + 773.5  cm’1 

6x  10" 18  0/300)  !/2 

6 

N2(l)  +  0  »  N2(0)  +  O  +  2329.3  cm'1 

3.2  x  10* 15  O/300)2-6 

7 

N2(0)  +  O(!D)  <=>  N2(v)  +  0(3P) 

2.4  x  10"11 

8 

>C02(10012)  +  N2(0)  <=>  'C02(10002)  +  N2(l)  -  2.5  cm'1 

5.0  x  10"13(300/T)1/2 

9 

•C02(02211)  +  N2(0)  <=>  1C02(02201)  +  N2(l)  -  5.8  cm"1 

5.0  x  10"13(300/T)1/2 

10 

1  C02(  1 001 1)  +  N2(0)  »  1C02(10001)  +  N2(l)  -  3.3  cm'1 

5.0  x  10"13(300/T)1/2 

11 

/CO2(01  111)  +  N2(0)  <=>  /CO2(01 101 )  +  N2(l)  +  6.7  cm’1 

5.0  x  1 0" 1 3(300/T) 1/2 

12 

'CO2(01111)  + M  <=>  'C02(04401)  +  M  +  332.3  cm'1 

2.22  x  lO'15 

+  1.22  x  10- 10  exp(-76.75/T°-33) 

13 

'CO2(0 1111)  +  O2(0)  <=>  ,C02(0220 1 )  +  02(  1 )  +  1 1 2.5  cm" 1 

3  x  10‘15  (1  +0.020  -  210)) 

14 

'CO2(01111)  +  O  <=>  'C02(04401)  + 0  +  332.3  cm"1 

2  x  lO' 13  0/300) >/2 

a  The  superscript  i  refers  to  each  of  the  isotopic  designations  626,  636,  628,  and  627.  In  reactions  (2)  and  (12), 
M  refers  to  N2  and  O2.  The  exothermicity  stated  is  for  reactions  involving  the  major  isotope  only. 

b  The  rate'' constants  are  for  the  reactions  in  the  forward  direction  and  are  the  same  as  those  used  by  Lopez-Puertas 
al.  [1986b].  T  is  the  kinetic  temperature  in  kelvin. 


02201  and  03301  states  are  similar  to  those  shown  in  Figure  6 
of  Wintersieiner  et  al.  [1992]. 

In  separate  steps,  calculations  are  then  performed  for  the 
higher-lying  states.  The  01111  populations  are  found  in  a 
calculation  that  uses  the  C02(V|V2)  populations  and  the 
previously  determined  N2(v=l)  population.  This  procedure  is 
valid  because  (1)  the  V]V2  manifold  populations  are  quite 
independent  of  v3  excitation,  as  discussed  in  the  previous 
paragraph,  and  (2)  the  N2(v=l)  population  is  determined 
primarily  by  exchange  of  vibrational  quanta  with  C02  states 
other  than  01111.  The  radiative  transfer  in  the  01101-01111 
band  is  calculated  using  the  RAD  code,  considering  the 
collisional  processes  11-14  of  Table  2.  This  has  been  done 
for  the  626  and  636  isotopes,  the  others  being  expected  to 
contribute  little  to  the  4.3-jim  radiance  in  the  hot  band. 

The  method  used  to  calculate  the  populations  of  the  higher- 
lying  states  of  Figure  1  (above  3500  cm'1,  groups  1,  2,  and  3) 
under  daytime  conditions  has  been  described  previously 
[Sharma  and  Wintersteiner ,  1985].  These  states  are  very 
effectively  populated  in  the  daytime  by  the  intense  short-wave 
solar  flux  at  2.7  and  2.0  pm,  but  relatively  little  of  the  energy 
is  reemitted  at  these  wavelengths  because  of  the  strong 
competition  from  the  4.3-jJ.m  emission  and,  at  lower  altitudes, 
the  V-V  exchange  with  N2  (e.g.,  reactions  (8)-(10)  of  Table  2). 
Thus  the  overall  short-wave  flux  is  almost  entirely  solar  in 
origin,  as  discussed  by  Sharma  and  Wintersteiner  [1985]. 
Also  described  in  that  paper  are  the  collisional  interactions 
among  the  upper  C02  states  of  reactions  (8)-(10);  these 
interactions  are  included  in  the  determination  of  the 


populations  of  those  states.  A  similar  statement  may  be  made 
for  the  group  2  and  group  3  states.  Because  of  the  dominance 
of  solar  pumping  in  the  excitation  of  all  of  these  higher-lying 
states,  the  calculation  neglects  the  airglow  (earthshine) 
absorption  contributions  to  the  populations  of  these  states  in 
the  daytime. 

2.2.  Model  Atmosphere  for  the  Spectral  Infrared 
Rocket  Experiment  (SPIRE)  Simulation 

To  calculate  the  populations  of  the  4.3-pm  states,  a  model 
atmosphere  specifying  temperature,  [C02],  [N2],  [02],  [O],  and 
[0(lD)]  as  functions  of  altitude  is  required.  The  basic  model 
atmosphere  used  here  is  the  same  as  that  described  by 
Wintersteiner  et  al.  [1992]  for  application  to  the  I5-pm 
bands.  For  temperature  and  the  major  constituents  the  model 
relies  heavily  on  MSIS  86  [Hedin,  1987]  above  the  mesopause 
and  the  FASCODE  midwinter  subarctic  atmosphere  [Anderson 
et  al.,  1986]  below  the  mesopause.  The  CCb  mixing  ratio 
profile  used  here  is  derived  from  experimental  data  and  is 
described  as  model  D  by  Wintersteiner  et  al.  [1992].  Atomic 
oxygen  and  C02  number  densities  for  this  model  are  shown  in 
Figure  2  of  that  paper.  The  0(1D)  number  densities  are  taken 
from  a  one-dimensional  diurnal  photochemical  model  which 
was  run  for  the  conditions  of  the  sunlit  SPIRE  data  (U. 
Makhiouf,  private  communication,  1993). 

2.3.  SPIRE  Data  Set 

The  Spectral  Infrared  Rocket  Experiment  (SPIRE)  observed 
the  earthlimb  in  the  1.4-  to  16.5-|im  spectral  region  with  a 
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circular  variable  filter  (CVF)  spectrometer  from  the  vicinity  of 
Poker  Flat,  Alaska,  at  dawn  [Stair  et  al.,  1985].  There  were  12 
spatial  limb  scans,  each  consisting  of  a  number  of  spectral 
scans,  with  tangent  heights  ranging  from  0  to  200  km.  The 
payload  was  exoatmospheric  and  sunlit.  Seven  spatial  scans 
looked  across  the  dawn  terminator,  the  eighth  looked  almost 
entirely  at  darkened  regions,  and  for  the  other  four  spatial 
scans  the  lines-of-sight  were  entirely  sunlit.  These  three 
groups  of  scans  are  herein  referred  to  as  terminator,  night,  and 
day  scans,  respectively.  The  principal  results  of  the 
experiment  were  summarized  by  Stair  et  ai  [1985]. 

The  experimental  4.3-|im  radiance,  integrated  from  4.12  to 
4.49  jim,  is  shown  in  Figure  2.  Scans  1-8  are  the  terminator 
and  night  scans  referred  to  above,  where  the  instrument  looked 
generally  toward  the  west  or  northwest.  Scans  9-12  are  the 
sunlit  day  scans  where  the  instrument  looked  generally  toward 
the  northeast  and  east  [Stair  et  al.,  1985].  The  great  disparity 
between  these  two  data  sets  emphasizes  the  important  role  that 
solar  pumping  plays  in  exciting  the  C02  states  that  emit  at 
4.3  pm.  For  the  sunlit  scans,  the  lines-of-sight  were 
illuminated  by  the  Sun  over  their  full  extent.  Nevertheless, 
the  scans  were  essentially  for  "low-Sun"  conditions.  The  solar 
zenith  angle  (SZA)  varied  along  each  viewing  path  and  varied 
in  a  different  way  for  different  spatial  scans.  For  those 
spectral  scans  with  tangent  heights  near  80  km  the  SZA  at  the 
tangent  point  varied  from  78°  to  84°  for  spatial  scans  9-12. 
Similar  variations  in  SZA  occur  for  other  tangent  heights  (e.g. 
77°  to  83°  for  40  km  and  79°  to  86°  for  120  km).  There  is  no 
apparent  systematic  variation  in  the  observed  limb  radiance  in 
different  spatial  scans  9-12  for  a  given  tangent  height  (see 
Figure  1 1).  .These  day  scans  are  therefore  modeled  with  a  fixed 
,tSZA  of  78°. 

3.  Results  for  the  SPIRE  Simulation 

3.1.  Vibrational  Temperatures 

Vibrational  temperatures  are  presented  for  both  nighttime 
(no  solar  excitation)  and  low-Sun  daytime  conditions  to 
emphasize  the  importance  of  solar  excitation  processes.  In 
this  work,  all  vibrational  temperatures  are  defined  with  respect 
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Figure  2.  Spectral  Infrared  Rocket  Experiment  (SPIRE) 
Hmb  radiance  in  the  band  4.12-4.49  jim.  as  a  function  of 
tangent  height.  Terminator  and  night  scans  1-8  are 
distinguished  from  sunlit  day  scans  9-12. 


to  the  ground  state,  as  indicated  in  equation  (1)  of 
Wintersteiner  et  al .  [1992].  Figure  3  shows  the  calculated 
vibrational  temperatures  of  the  C02(00011)  state  for  the  four 
isotopes  studied  here,  using  the  model  atmosphere  discussed 
above  and  assuming  no  solar  excitation.  Figure  4  gives  the 
same  vibrational  temperatures  for  low-Sun  daytime  conditions 
with  a  solar  zenith  angle  of  78°.  Unlike  the  C02(v2)  case, 
where  the  main  (626)  isotope  may  be  close  to  LTE  as  high  as 
100  km  [Lopez-Puertas  et  al. ,  1992],  the  v3  states  deviate  from 
LTE  near  60  km  for  nighttime  conditions  and  near  45  km  for 
daytime  (sunshine-pumped)  conditions.  From  60  km  to  near 
100  km  (night)  and  from  45  km  to  near  110  km  (day)  the 
vibrational  temperature  Tvjb  is  greater  than  the  kinetic 
temperature  T,  as  radiative  excitation  from  sunshine  or  from 
upwelling  earthshine  from  the  warm  stratosphere  becomes  an 
important  excitation  process.  Above  100-110  km.  Tvib  falls 
below  T  because  Tvjb  tends  to  a  radiative  equilibrium  profile 
and  is  thus  independent  of  the  local  kinetic  temperature. 

Solar  flux  absorption  in  the  4.3-|im  bands  accounts  for  the 
large  differences  between  the  daytime  results  of  Figure  4  and 
the  nighttime  results  shown  in  Figure  3.  It  is  responsible  for 
most  features  of  the  daytime  profiles  above  60  km.  In 
addition,  the  impact  of  O (lD)  on  the  C02  vibrational 
temperatures  is  greater  than  previously  reported.  Without  the 
quenching  of  0(lD)  by  N2  the  vibrational  temperatures  shown 
in  Figure  4  would  be  lower  throughout  the  mesosphere  and 
lower  thermosphere,  and  the  curves  would  not  deviate  from 
LTE  until  55-60  km.  This  will  be  discussed  further  in  section 
4.1.  The  differences  between  the  isotopic  vibrational 
temperatures  in  the  mesosphere  are  largely  due  to  the  fact  that 
the  minor  isotope  bands  have  smaller  opacity  than  the  band  of 
the  principal  (626)  isotope  due  to  much  smaller  number 
densities.  This  means  that  the  solar  flux  penetrates  to  lower 
altitudes  in  the  minor  species  bands  in  the  daytime  and  the 
stratospheric  earthshine  penetrates  more  effectively  into  the 
mesopause  region  at  night.  This  results  in  higher  values  of 
Tvjb  for  the  minor  isotope  levels.  Below  100  km  the  00011 
vibrational  temperatures  can  be  much  higher  than  are  shown  in 
Figure  4  if  the  solar  zenith  angle  (SZA)  is  smaller.  Because  of 
radiative  transfer  upward  from  this  region,  the  vibrational 
temperatures  of  the  00011  states  at  higher  altitudes  (110  km 
and  above)  may  thus  depend  on  the  SZA  even  though  the  direct 
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Figure  3.  Nighttime  vibrational  temperatures  for  the 
00011  state  for  four  C02  isotopes.  The  kinetic  temperature  is 
also  shown  for  purposes  of  comparison. 
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Figure  4.  Vibrational  temperatures  for  the  0001 1  state  for 
four  C02  isotopes,  for  daytime  conditions  with  a  solar  zenith 
angle  of  78°.  The  kinetic  temperature  is  also  shown. 


solar  excitation  rate  does  not.  This  point  will  be  discussed 
further  in  section  4.1. 

Figure  5  shows  the  N2  vibrational  temperatures  from  the 
same  calculations  that  produced  Figures  3  and  4.  along  with 
the  corresponding  626  C02(00011)  vibrational  temperatures. 
At  the  lower  altitudes  the  N2  and  CCA  states  are  very  closely 
coupled  and  the  N2  vibrational  temperature  approaches  the 
C02  vibrational  temperature  for  both  day  and  night  conditions 
because  of  the  predominance  of  the  rapid  intermolecular  V-V 
process;  equation  (1).  At  the  higher  altitudes,  on  the  other 
hand,  the  N2  vibrational  temperature  approaches  the  kinetic 
temperature  at  night,  due  to  the  importance  of  the  thermal 
process  whereby  oxygen  atoms  directly  excite  N2  vibration 
(Table  2,  reaction  (6)).  The  N2(v=l)  lifetime  is  determined  in 
this  region  by  collisions  with  atomic  oxygen  rather  than  the 
V-V  process,  equation  (1);  thus  the  N2  vibrational  temperature 
decouples  from  .the  C02(00011)  vibrational  temperature  and 
couples  instead  to  the  kinetic  temperature  at  night.  Under 
sunlit  conditions  the  N2  vibrational  temperature  also 
decouples  from  the  C02(00011)  vibrational  temperature  at  the 
higher  altitudes  but  in  this  case  is  driven  considerably  above 
the  kinetic  temperature  because  of  the  direct  quenching  of 
O(’D)  by  N2. 

The  vibrational  temperature  profiles  for  the  CO2(0 1111) 
state  for  the  626  and  636  isotopes,  for  daytime  conditions 
with  a  solar  zenith  angle  of  78°,  are  shown  in  Figure  6.  The 
shape  of  the  curves  in  the  lower  mesosphere  is  similar  to  that 
of  N2(v=1),  indicating  the  influence  of  V-V  transfer  as  an 
excitation  process.  At  70  km  and  higher,  radiative  processes 
are  the  principal  factors  determining  the  01111  temperatures. 
The  differences  between  the  626  and  636  curves  are  largely  due 
to  the  different  optical  thicknesses  of  the  respective  bands: 
for  the  636  isotope  there  is  deeper  penetration  of  solar  flux, 
and  a  greater  fraction  of  the  upwelling  radiation  reaches  the 
middle  mesosphere.  In  the  thermosphere  the  vibrational 
temperatures  of  this  level  are  larger  than  those  of  previous 
calculations  (see,  for  example,  Lopez-Puertas  and  Taylor 
[1989],  Figure  5).  This  is  due  to  increased  population  of  the 
01101  level  (compared  to  previous  results)  which  in  turn  is 
caused  by  the  greater  rate  constant  for  collisions  between 
C02(v2)  and  atomic  oxygen  [Wintersteiner  et  al. ,  1992]. 


Figures  7-10  show  the  vibrational  temperatures  of  the 
higher-lying  groups  1,  2,  and  3  states  of  Figure  1  (-3600- 
5000  cirf1).  The  values  of  Tvib  for  the  02211  and  03311 
states  fall  progressively  further  below  Tvib  values  for  the  other 
states  in  their  respective  groups  as  altitude  increases  above  60 
km.  This  is  because  they  are  not  effectively  pumped  by  solar 
flux,  and  the  collisional  coupling  with  the  other  states  which 
are  pumped  at  2.7  |im  becomes  inefficient  in  the  mesosphere 
and  thermosphere. 

3.2.  Limb  Radiance 

The  daytime  scans  of  the  SPIRE  experiment  were  modeled 
by  using  the  vibrational  temperatures  shown  in  Figures  4  and 
6-10,  calculated  for  a  solar  zenith  angle  of  78°,  as  input  to  the 
line-by-line  radiance  code  NLTE  [Wintersteiner  and  Sharma , 
1985].  This  code  calculates  integrated  radiance  in  a  limb  view 
for  each  spectral  line  within  a  band  for  selected  tangent 
heights  and  then  sums  the  results  to  obtain  total  band 
radiance.  The  calculation  must  be  done  for  each  band  and  each 
isotope  separately.  The  results,  using  all  the  bands  listed  in 
Table  1,  are  given  in  Figure  1 1,  along  with  the  data  points  for 
the  day  lit  scans  9-12  which  are  repeated  from  Figure  2.  The 
nominal  instrument  noise  level  of  about  4xl0"9  Watt  cm"2  sr’1 
for  this  wavelength  region  [S/a/r  <>/  aL ,  1985]  has  been 
subtracted  from  the  data  values  for  this  figure.  The  agreement 
of  our  model  with  the  data  at  all  altitudes  is  regarded  as 
validation  of  the  model  and  confirmation  that  it  properly 
accounts  for  the  processes  that  are  important  for  exciting 
C02(v3)  states  in  the  daytime. 

Figure  12  gives  the  breakdown  of  the  daytime  model 
radiance  into  contributions  from  individual  bands  or  groups  of 
bands.  Several  features  are  noteworthy:  (1)  The  strong  626 
fundamental  00011-00001  band  is  responsible  for  virtually  all 
of  the  radiance  above  1 10  km.  (2)  The  626  fundamental  plays 
a  relatively  unimportant  role  for  limb  paths  through  the 
mesosphere  (55-85  km)  where  the  strong  band  has  become 
self-absorbed.  (3)  The  major  contribution  to  the  total  limb 
radiance  for  tangent  heights  of  55-95  km  is  made  by  the  group 
1  fluorescent  states  lying  near  3500  cm'1  (see  Figure  1), 
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Figure  5.  N2  vibrational  temperatures  for  nighttime 
conditions  and  for  daytime  conditions  with  a  solar  zenith 
angle  of  78  .  The  kinetic  temperature  is  shown  for  purposes 
of  comparison,  as  are  the  night  and  day  C02(00011) 
vibrational  temperatures  for  the  626  isotope. 
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Figure  6.  CO2(0 1111)  vibrational  temperatures  for 

daytime  conditions  with  a  solar  zenith  angle  of  78°.  The 
kinetic  temperature  is  included  for  purposes  of  comparison. 
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Figure  8.  Vibrational  temperatures  of  the  group  2  states 
11112,  11111,  and  03311  of  the  626  isotope  for  daytime 
conditions  with  a  solar  zenith  angle  of  78°. 


particularly  the  10011-10001  and  10012-10002  bands,  which 
for  the  80-km  tangent  path  account  for  60%  of  the  observed 
limb  radiance.  Below  80  km  the  radiance  contribution  from 
the  group  1  bands  falls  off  as  they  become  self-absorbed.  (4) 
For  lower  tangent  heights  near  50  km  the  contributions  from 
all  the  groups  of  bands  are  roughly  comparable,  reflecting  the 
differing  optical  thicknesses  of  the  various  bands.  The 
dominance  of  the  fluorescent  states  from  55  to  95  km  and  the 
relative  unimportance  of  the  626  fundamental  compared  even 
to  the  minor  isotope  fundamentals  in  the  mesosphere  were 
'I  predicted  by  Kumer  [19771  based  on  a  band  model.  Similar 
"conclusions  were  reached  by  Lopez-Puenas  and  Taylor  [1989] 
from  SAMS  measurements.  The  present  work  has  reached 
these  conclusions  using  a  full  line-by-line  treatment  and  is 
validated  by  good  agreement  with  the  SPIRE  data  set. 

4.  Discussion 

4.1.  Comparison  With  Previous  Studies 

The  limb  radiance  profiles  obtained  here  for  the  different 
bands  (see  Figure  12)  are  similar  in  general  appearance  to 
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Figure  7.  Vibrational  temperatures  of  the  group  1  states 
10012,  10011,  and  02211  of  the  626  isotope  for  daytime 
conditions  with  a  solar  zenith  angle  of  78°. 


those  calculated  by  Kumer  and  James  [1982],  but  there  are 
differences  in  the  details.  The  present  calculations  predict  a 
lower  limb  radiance  in  the  vicinity  of  80  km  for  scans  having 
comparable  solar  elevation  and  a  more  rapid  drop-off  with 
increasing  tangent  altitude  above  80  km.  Some  of  this  results 
from  substantially  smaller  contributions  from  the  01  111 
bands  and  the  minor  isotope  fundamentals  (collectively,  the 
"weak”  bands  in  their  report)  which  occur  throughout  the  range 
of  tangent  heights  and  is  most  likely  due  to  the  great 
improvement  in  the  treatment  of  radiative  transfer  in  these 
bands.  The  present  model  also  predicts  smaller  contributions 
from  the  626  fundamental  for  tangent  heights  below  75  km, 
which  is  probably  because  the  line-by-line  approach  enables  a 
very  careful  treatment  of  the  Voigt  wings  of  this  extremely 
opaque  band.  However,  we  include  bands  originating  in  group 
2  and  group  3  which  contribute  below  75  km;  these  are  not  in 
the  earlier  model,  and  they  partially  compensate  for  the 
discrepancies  in  the  other  bands  noted  above. 

Lopez-Puertas  and  Taylor  [1989],  in  analyzing  SAMS  data, 
also  found  that  emission  in  hot  bands  dominates  the  limb 
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Figure  9.  Vibrational  temperatures  of  the  group  3  states 
20013.  20012.  20011,  12212,  12211,  and  04411  of  the  626 
isotope  for  daytime  conditions  with  a  solar  zenith  angle  of 
78°. 


NEBEL  ET  AL.:  CCb  4.3  |im  DAYGLOW 


10.416 

✓ 


GROUP  I  (636) 


VIBRATIONAL  TEMPERATURE  (K) 


Figure  10.  Vibrational  temperatures  of  the  group  1  states 
as  in  Figure  7,  but  for  the  636  isotope. 

radiance  in  the  mesopause  region,  but  a  direct  comparison  with 
the  present  study  is  frustrated  by  the  fact  that  the  SAMS 
instrument  uses  the  pressure  modulation  technique  which 
makes  it  much  more  sensitive  to  the  stronger  bands.  In 
addition,  for  the  70°-50°N  latitude  range,  they  include  the 
01111-01101  bands  and  the  group  1- bands  in  a  single  profile, 
making  the  determination  of  relative  contributions  from  these 
bands  impossible.  For  lower  latitudes  they  do  find  that  the 
group  1  bands  give  a  larger  contribution  than  the  01111- 
01101  bands  inythe  mesopause  region,  but  the  fundamental 
bands  are  also  important  in  this  region  in  their  study.  They 
al&o^use  a  different  temperature  profile  which  further 
complicates  a  direct  comparison  with  the  present  study. 

As  mentioned  earlier,  the  quenching  of  0(lD)  by  N2  has  a 
greater  impact  on  the  C02  vibrational  temperatures  than  that 
reported  previously  by  other  authors.  Kumer  [1977] 
considered  this  process  and  obtained  a  maximum  enhancement 
of  the  C02(00011)  vibrational  temperature  of  ~4  K  at  62  km, 
using  a  solar  zenith  angle  of  60°.  We  find  an  enhancement  of 
~5  K  for  the  major  626  isotope  due  to  this  process  at  the  same 
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Figure  11.  Comparison  of  the  SPIRE  limb  radiance  in  the 
band  4.12-4.49  pm  with  the  daytime  model  calculation.  Data 
from  sunlit  scans  9-12  only  are  plotted. 
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Figure  12.  Breakdown  of  the  daytime  model  radiance  into 
contributions  from  individual  bands  or  groups  of  bands. 
Radiance  from  the  01111  band  and  the  group  1  bands  include 
contributions  from  the  626  and  636  isotopes.  Radiance  from 
the  group  2  and  group  3  bands  include  contributions  from  the 
626  isotope  only. 


altitude,  at  a  solar  zenith  angle  of  78°.  Harris  and  Adams 
[1983]  found  that  the  CCM00011)  concentration  is  essentially 
unaffected  by  the  addition  of  0(lD),  and  Lopcz-Puerias  et  ai 
[1986b]  reported  that  the  C02(00011)  population  is  not 
significantly  perturbed  by  the  quenching  of  0(]D)  by  N2.  at 
least  for  overhead  Sun  conditions.  In  contrast,  we  find  that  the 
C02(0001 1)  vibrational  temperature  is  enhanced  by  4-6  K  due 
to  this  process  throughout  the  mesosphere,  with  comparable 
enhancements  in  the  lower  thermosphere. 

As  indicated  in  Figure  5.  the  vibrational  temperature  of  the 
626  fundamental  is  about  298  K  in  the  thermosphere  above 
110  km.  Previous  calculations  have  yielded  values  between 
300  and  310  K  [Kumer,  1977;  Shved  et  ai,  1978:  Lopez - 
Puertas  et  a!..  1986b;  Lopez-Puertas  and  Taylor.  1989].  The 
reason  for  this  difference  is  that  all  of  the  previous 
calculations  were  done  for  a  higher  Sun  than  is  used  in  the 
present  set  of  calculations.  When  the  Sun  is  higher,  the 
CC>2(00011)  state  will  be  populated  above  110  km  not  only  by 
direct  solar  excitation  but  also  by  an  increased  amount  of 
upwelling  radiation  (earthshine)  from  lower  altitudes.  The 
*  latter  is  due  to  emission  from  molecules  excited  by  sunlight 
that  penetrates  more  deeply  into  the  atmosphere  than  it  does 
for  the  low  Sun  angle  considered  in  the  present  study.  Thus 
earthshine  from  below  is  relatively  more  important  for  a  high 
Sun  than  it  is  for  a  low  Sun  in  determining  the  vibrational 
temperature  of  the  C02(00011)  state  in  the  thermosphere. 

This  is  indicated  in  Figure  13  which  shows  the  626(00011) 
and  N2(v=l)  vibrational  temperatures  for  solar  zenith  angles 
of  0°,  55°,  and  78°.  These  curves  have  been  generated  using  a 
single  profile  for  0(lD)  (appropriate  for  a  solar  zenith  angle  of 
78°)  because  the  purpose  of  the  curves  is  to  compare  the 
importance  of  upwelling  radiation  from  below  on  the 
C02(00011)  vibrational  temperatures  in  the  thermosphere  for 
different  solar  zenith  angles.  Note  the  large  variation  in  the 
CO 2  curves  in  the  80-  to  110-km  region  due  to  the  deeper 
penetration  of  sunlight  for  the  lower  solar  zenith  angles  ! 

(higher  Sun). 
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N2  and  C02  VIBRATIONAL  TEMPERATURES 


Figure  13.  Vibrational  temperatures  of  the  N2(v=l)  state 
and  the  00011  state  of  C02  (626  isotope)  for  solar  zenith 
angles  of  0°,  55°,  and  78°. 


4.2.  Approximations  and  Uncertainties 

Several  approximations  and  uncertainties  which  are 
inherent  in  the  present  set  of  calculations  are  now  discussed. 
First  of  all.  there  are  factors  related  to  the  precision  of  the 
calculations  themselves.  Questions  of  convergence,  line 
overlap,  layering,  and  choice  of  the  lower  boundary  for  the 
RAD  algorithm  have  been  discussed  previously  [Wintersteiner 
ff  et  dl.,  1992].  As  argued  there,  any  inaccuracies  resulting  from 
these  factors  would  be  very  small  compared  with  some  of  the 
uncertainties  discussed  below. 

It  has  been  assumed  here  that  there  are  no  day-night 
differences  in  the  populations  of  the  02201  and  03301  states 
[Wintersteiner  et  al. .  1992].  Lopez-Puertas  et  al.  [1992]  have 
found  that  the  day-night  differences  in  the  vibrational 
temperatures  of  these  states  may  be  as  much  as  8  K  for  some 
atmospheric  conditions  for  the  02201  state  and  14  K  for  the 
03301  state  in  the  80-  to  90-km  altitude  range.  To  test 
whether  the  increased  opacity  resulting  from  larger  lower-state 
populations  would  diminish  the  limb  radiance  obtained  for  the 
group  1  bands,  a  set  of  ad  hoc  2v2  vibrational  temperature 
profiles  was  created  with  an  8  K  enhancement  near  the 
mesopause.  and  the  limb  radiance  was  recalculated.  The  result 
is  a  1%  decrease  in  radiance  in  the  group  1  bands  at  85  km. 
with  smaller  differences  for  other  paths.  The  group  2  bands  are 
optically  thin,  so  no  change  is  possible  due  to  higher  3v-> 
populations.  Thus  the  effect  of  possible  dav-night  differences 
in  the  populations  of  states  02201  and  03301  is  negligible  for 
the  present  study. 

Lopez-Puertas  et  al.  [1986b]  (hereinafter  referred  to  as  LPb) 
consider  a  number  of  coliisional  processes  which  are  not 
included  here.  These  processes  include  isotope  exchange  of  v3 
quanta  (reaction  (2)  from  Table  2  of  LPb)  and  the  coliisional 
relaxation  processes  of  C02  (0001 1)  on  N2,  02,  and  O  leading 
to  the  C02  (02201)  state  (reactions  (3b),  (4b),  and  (7b)  from 
the  same  table).  The  primary  justification  for  omitting  these 
processes  is  that  the  present  work  treats  daytime  limb  radiance 
only,  where  solar  pumping  is  clearly  the  dominant  excitation 
process  for  the  4.3-Jim  emitting  states.  Additional  processes 
considered  by  LPb  are  expected  to  change  the  4.3-|im  radiance 
by  very  little  for  the  sunlit  case,  and  in  the  interest  of 
simplicity  they  are  neglected  in  this  work. 


In  addition,  LPb  consider  excitation  of  N2(v=l)  by  OH*(v  < 

9)  (reaction  (9)  from  Table  2  of  LPb).  OH,  which  is  produced 
in  excited  vibrational  states  (v  <  9)  near  the  mesopause  by  the 
reaction  of  H  and  O3,  may  affect  the  4.3-|im  radiance  in  two 
ways:  (1)  by  exchanging  energy  with  the  N2-C02  reservoir 
through  V-V  transfer  [Kumer  et  al .,  1978]  and  (2)  by  direct 
emission  in  the  9-8  and  8-7  bands.  Since  the  OH  production 
rate  is  small  in  the  daytime  due  to  the  photolysis  of  03,  these 
contributions  have  been  neglected  in  the  calculations  reported 
above.  Even  at  night  it  is  found  that  the  enhancements  of  the 
inband  4.3-|im  limb  radiance  due  to  OH  are  only  -2-5  x  10-9 
Watt  cm"2  sr"1  at  most,  so  this  omission  is  of  no  consequence 
here. 

Reactions  (12)-(16)  from  Table  2  of  LPb  involve  the  state 
02(v=l),  the  population  of  which  is  considered  an  unknown 
quantity  in  their  investigation.  It  is  assumed  here  that  the 
02(v=l)  state  is  in  LTE.  This  is  certainly  the  case  above  85 
km  where  the  V-T  process  with  atomic  oxygen  dominates  the 
behavior  of  02(v=l);  thus  the  02(v=l)  vibrational  temperature 
is  coupled  to  the  kinetic  temperature  T.  Below  45  km, 
everything  is  in  LTE,  including  V-V  processes,  since 
collisions  are  so  frequent.  Thus  it  is  only  in  the  mesosphere 
(50-  to  80-km  altitude  range)  that  02(v=l)  may  not  be  in  LTE. 
We  did  a  simple  perturbation  calculation  of  the  02(v=l) 
vibrational  temperature,  holding  other  quantities  Fixed,  and 
determined  that  it  exceeds  the  local  kinetic  temperature  by 
amounts  ranging  from  3.6  K  to  less  than  1  K  over  this  altitude 
range.  We  then  calculated  the  additional  excitation  of  N2  that 
would  occur  through  reaction  (5)  (see  Table  2)  if  these  upper 
bound  estimates  were  correct  and  found  that  it  is  between  20 
and  200  times  smaller  than  that  due  to  0(!D).  The  relaxation 
of  the  LTE  restriction  on  the  02(v=l)  vibrational  temperature 
could  therefore  enhance  the  C02(00011)  vibrational 
temperatures,  which  are  strongly  coupled  to  N2,  by  only  a 
small  fraction  of  a  degree  in  the  mesosphere.  Thus  the 
resulting  limb  radiance  is  essentially  unaffected  by  the 
assumption  of  LTE  for  the  02(v=l)  state. 

As  pointed  out  by  Wintersteiner  et  ai  [1992],  there  is 
considerable  variability  in  the  C02  mixing  ratio  profile  based 
on  experimental  information  from  rocket-  and  satellite-borne 
instruments  (see  Figure  3  of  that  paper).  Although  the 
vibrational  temperatures  calculated  here  are  quite  insensitive 
to  the  C02  densities  [Wintersteiner  et  al. .  1990],  these 
densities  will  have  an  appreciable  effect  on  total  limb-view 
radiance,  especially  at  the  higher  altitudes  where  the  limb 
radiance  is  directly  proportional  to  C02  column  density  along 
the  viewing  path.  The  C02  mixing  ratio  profile  used  here 
represents  an  average  of  data  reported  by  rocket  experiments 
in  the  literature  (curve  D  of  Figure  2  of  Wintersteiner  et  al. 
[1992]).  Error  bars  on  the  experimental  mixing  ratio  values 
are  of  the  order  of  a  factor  of  2  at  the  higher  altitudes:  thus  the 
slight  underprediction  of  radiance  values  for  110-130  km  seen 
in  Figure  1 1  is  not  significant. 

5.  Summary  and  Conclusions 

The  RAD  code  calculates  excited  state  population  densities 
under  non-LTE  conditions  in  the  atmosphere  using  line-by¬ 
line  radiative  transfer  algorithms.  The  line-by-line 
algorithms  have  been  used  to  calculate  population  densities 
and  the  corresponding  vibrational  temperature  profiles  for  the 
important  4.3-pm  emitting  states  of  CO->  under  daytime 
conditions.  The  quenching  of  0(‘D)  by  N2  has  a  greater  impact 


on  these  vibrational  temperatures  than  previously  reported  by 
other  authors.  Total  integrated  radiance  in  a  limb  view  from 
the  bands  associated  with  each  of  these  states  has  also  been 
calculated  using  the  NLTE  code.  Solar  pumping  is  the 
dominant  excitation  process  for  the  4.3-jim  emitting  states  in 
the  daytime:  the  total  limb  radiance  is  not  very  sensitive  to 
any  other  excitation  processes.  The  largest  contributors  to 
the  radiance  for  tangent  heights  in  the  upper  mesosphere  and 
lower  thermosphere  (55-95  km)  are  the  v3  bands  10012-10002 
and  10011-10001  (see  Figure  1)  where  the  upper  states  absorb 
sunlight  at  2.7  Jim  and  reemit  preferentially  at  4.3  fim.  The 
agreement  of  the  model  calculations  with  SPIRE  sunlit  data  at 
all  altitudes,  as  shown  in  Figure  11,  indicates  that  the  model 
properly  accounts  for  the  important  excitation  processes  of 
C02(v3)  in  the  daytime.  Thus  the  model  has  been  validated  by 
agreement  with  experimental  results,  indicating  that  the 
approach  is  sound.  Future  work  will  concentrate  on  the  night 
and  terminator  models,  on  comparison  of  the  model  with  other 
data  sets,  and  on  adding  other  processes  to  the  model. 
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